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Summary 


The basic objective of this research was to identify, develop and recommend 
turbulence models which could be incorporated into CFD codes used in the design of the 
NASP vehicles. To accomplish this goal, a combined effort consisting of experimental and 
theoretical phases was undertaken. The experimental phase consisted of a literature survey 
to collect and assess a database of well documented experimental flows, with emphasis on 
high speed or hypersonic flows, which could be used to validate turbulence models. Since 
it was anticipated that this database would be incomplete and would need supplementing, 
additional experiments in the NASA Ames 3.5' hypersonic wind tunnel (HWT) were also 
undertaken. The theoretical phase consisted of identifying promising turbulence models 
through applications to simple flows, and then investigating more promising models in 
applications to complex flows. The complex flows were selected from the database devel- 
oped in the first phase of the study. For these flows it was anticipated that model perfor- 
mance would not be entirely satisfactory so that model improvements or corrections 
would be required. The primary goals of the investigation were essentially achieved. A 
large database of flows was collected and assessed, a number of additional hypersonic 
experiments were conducted in the Ames HWT, and two turbulence models (£-e and k-ca 
models with corrections) were determined which gave superior performances for most of 
the flows studied and are now recommended for NASP applications. 

Introduction 

With the advent of the hypersonic airplane, hypersonic flows are receiving spe- 
cial attention from researchers in computational fluid dynamics. Complex flow phenom- 
ena such as shock-wave boundary layer interactions, separation and combustion are of 
particular interest because of their importance to the successful design of structural, pro- 
pulsive and thermal protection systems. Rapid advances in CFD in recent years have 
resulted in its increased use as a design tool for aeronautical systems and have lead to 
reductions in the time and costs of wind tunnel testing. A major obstacle to the use of CFD 
as a design tool is its dependence on turbulence modeling for accurate prediction of com- 
plex flows. Although recent advances have been made in turbulence modeling, many more 
will be required before CFD can be applied with confidence to a wide range of flow prob- 
lems. This is especially true at hypersonic speeds where high temperatures and pressures 
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create additional difficulties for turbulence modeling and where wind tunnel experiments 
which can be used to validate models are sparse. 

The overall goal of this research is to identify, develop and recommend turbu- 
lence models which can be incorporated into the CFD codes used in the design of the 
NASP vehicle. To accomplish this objective, a combined and joint effort consisting of 
experimental and theoretical phases was undertaken. The objective of the experimental 
phase was to conduct a literature survey to identify and assess a database of well docu- 
mented experimental flows which could be used to validate turbulence models. Since it 
was anticipated that this database would be incomplete, it was also decided to perform 
additional experiments in the NASA Ames 3.5’ hypersonic wind tunnel (HWT) to supple- 
ment the database. 

The objective of the theoretical phase was to identify promising turbulence mod- 
els through applications to simple flows, such as flat plate flows, and to then investigate 
the more promising models in applications to more complex flows. The complex flows 
were to be selected from the database developed in the first phase of the study. For these 
flows it was anticipated that model performance would not be entirely satisfactory so that 
model improvements or corrections would be required. The flows of interest were 
restricted to ideal gas flows because of the sparsity of high quality experimental validation 
data and viable turbulence models for real gas flows. 

The schedule of tasks and milestones for the completion of the research on Gov- 
ernment Work Package 18 is shown in fig. 1. With this report the work is essentially com- 
plete with the exception of the compressible shear layer experiment in the Ames 3.5’ 
HWT. Completion of this experiment was halted due to the lack of funds caused by fund- 
ing reductions in the NASP project. 

The report is organized into 7 sections. Following the introduction, two sections 
on the experimental phase of the study will be presented including one on the database 
collection and assessment and another on the experiments conducted in the Ames 3.5’ 
HWT. Next, sections describing the theoretical phase of the study will be presented. These 
include sections on recommended baseline turbulence models, compressibility corrections 
recommended for improved predictions of complex flows and representative results of 
numerical predictions using the baseline and corrected models. Finally, the report con- 
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eludes with a summary of basic results and recommendations including topics for future 
study. 


Database Collection and Assessment 

As stated in the introduction, the purpose of the database collection and assess- 
ment activity was to provide a base of reliable and well documented wind tunnel experi- 
ments by means of which turbulence models could be validated. These experiments, for 
the most part, involve relatively simple geometric shapes which may be viewed as sepa- 
rate elements of an overall vehicle. The results of this activity are reported by Settles and 
Dodson (1991, 1993a and 1993b) and will be summarized here. 

Settles and Dodson (1991 and 1993b) provide a survey and assessment of exper- 
iments involving two and three dimensional shock- wave boundary-layer interaction flows. 
Eight hundred experiments were initially identified for further review. Of these, 112 dis- 
tinct experiments were found involving flows at Mach 3 and above. The acceptance crite- 
ria applied to these included: 1) Measurements of surface pressures, skin friction and/or 
heat transfer, and velocity, temperature, or pitot pressure profiles at selected locations, 2) 
well defined experimental boundary conditions, 3) well defined error bounds, 4) adequate 
spatial resolution of measurements and 5) full documentation of tabulated data. For hyper- 
sonic conditions, i.e. M > 5, only 7 studies passed the acceptance criteria and only three 
were three dimensional. An additional 11 experiments passed the criteria at supersonic 
speeds. 


The survey and assessment of attached boundary layer and free shear flow exper- 
iments is given by Settles and Dodson (1993a). The acceptance criteria applied to these 
cases was identical to that applied to the shock-wave boundary-layer interaction experi- 
ments. For the boundary layer experiments, 153 candidate cases were identified for further 
review. Of these 39 were subjected to the acceptance criteria. No hypersonic and only 9 
supersonic cases passed. For the free shear layer experiments, 1137 candidate cases were 
identified for further review. Of these, 45 were subjected to the acceptance criteria and 
only 3 passed. 
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Experiments in the NASA Ames 3.5’ Hypersonic Wind 

Thnnel 

The experiments were conducted in the Ames 3.5’ hypersonic wind tunnel and 
were done with relatively simple generic shapes such as cones, cylinders, plates and 
wedges. These shapes were chosen to typify locations on a high speed vehicle where tur- 
bulence modeling was expected to be a critical issue. The experiments were run at nomi- 
nal Mach numbers ranging from 7 to 8.3 and unit Reynolds numbers (per meter) from 4.9 
to 5.8 million. Boundary-layers approaching the interaction region were relatively large, 
with thicknesses on the order of 2.5 to 3.7 cm, which allowed detailed flow field surveys to 
be easily made. In most of the experiments, both surface measurements and flow field 
(profile) measurements were made including initial boundary conditions required to start 
numerical computations. With two minor exceptions, the experiments passed the accep- 
tance criteria described in the previous section. In addition, an analysis of measurement 
errors was made and documented. The results of these experiments are described in detail 
by Kussoy et al (1989; 1991a, b; 1992; 1993a, b,c) and Horstman and Kussoy (1989) and 
the data have been made available on floppy disks. 

M = 7 Hypersonic Cylinder-Flare and Fin flows 

The test bed employed in this experiment consisted of an ogive-cylinder at zero 
angle of attack with a series of removable symmetric flares or sharp fins (see figs. (2a and 
2b)). Both flare and fin angles were varied, producing shock waves of various strengths, 
and resulting in both attached and separated flow fields. Detailed measurements verified a 
fully developed turbulent boundary- layer on the cylinder ahead of the interaction region. 
The resulting flows were axisymmetric with and without separation for the flare case, and 
three dimensional with separation for the fin case. Surface pressures and heat transfer rates 
were measured on both configurations, and flow field surveys were done on the flare con- 
figuration. The results are reported in Kussoy and Horstman (1989) and Horstman and 
Kussoy (1989). 


M = 8.2 2-D Wedge and 3-D Vertical Fin Flows 

A flat plate arrangement was used for this experimental series and is shown in 
fig.(3). It was of a hollow modular construction, enabling both test bodies and instrumen- 
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tation to be easily manipulated and changed. The full length of the wind tunnel test section 
was utilized in the design and because of this a well developed equilibrium turbulent 
boundary-layer was present at the shock interaction zone. Two configurations were tested; 
the first consisted of a sharp wedge supported over the width of the test section, the second 
was a sharp vertical fin attached to the plate surface. These are both illustrated in fig.(3). 
Both the wedge and fin angles were varied, producing shock-wave boundary-layer interac- 
tions of varying strength with a maximum wall-pressure ratio of p/poo = 21.5 and 6.4, 
respectively. This resulted in both attached and separated flow fields for the wedge flows, 
and swept three dimensional vortical flow fields for the fin flows. Detailed surveys verified 
a fully developed hypersonic turbulent boundary-layer on the flat plate ahead of the inter- 
action zone. 

For the wedge configuration, only surface conditions were measured. For the fin 
configuration, however, mean flow profile surveys were also taken - both in the undis- 
turbed and interaction regions - and from them pitot pressure contours and boundary-layer 
thickness parameters were obtained. We believe this is the first fully three dimensional 
shock- wave boundary-layer interaction flow to be so documented at hypersonic speeds. 
Results for both configurations are reported in Kussoy et al (1991a,b;1992). Experimental 
and computational results for the 10° and 15° vertical fin flows are discussed in the sec- 
tion on Model Validation-Representative Results. 

M = 8.3 Crossing Shock Flow 

For the third series of experiments, a configuration was chosen to reflect several 
key elements of a generic hypersonic inlet. These included a thick turbulent boundary- 
layer approaching two vertical fins of varying wedge angle, a crossing shock pattern pro- 
duced by the fins, boundary-layer vortices, large pressure gradients, and separation zones. 
The test body for this series of experiments is shown in fig.(4). 

Streamwise and transverse surface pressure and heat transfer distributions were 
measured as well as flow field surveys of pitot pressure and flow angle. One important 
result of these measurements should be mentioned here. This was the persistence of an 
extensive low pressure region far downstream of the fin leading edges. This low pressure 
region implied that the generic inlet tested here would not be a very efficient pressure dif- 
fusing device. The experimental results for this configuration are given in Kussoy et al 
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(1993a, b,c). experimental and computational results for the 15° fin-angle case are dis- 
cussed in the section on Model Validation-Representative Results. 

Recommended Baseline TXirbulence Models 

The baseline turbulence models recommended for use in NASP applications will 
be described in this section. A variety of turbulence models have been investigated in 
varying degrees throughout the course of the study. They include 0-eq, 1-eq, and 2-eq 
eddy viscosity models, and Reynolds stress transport models. For NASP applications, the 
primary emphasis has been placed on 2-eq models and these models will be the only ones 
described in detail here. Descriptions of other models used in the course of the study, and 
results obtained with them, are given in Coakley and Huang (1992) and Coakley and Mar- 
vin (1993), Horstman (1991,1992), and Huang and Coakley (1993a,b). 

The models investigated and recommended here are those that utilize no slip 
boundary conditions at solid walls and involve the use of wall damping functions (in most 
cases). This is in contrast to conventional practice with two equation models in which wall 
functions and essentially slip type boundary conditions are used. This choice was made 
because at the very high speeds of hypersonic flight, the effective Reynolds numbers of 
the flows can be quite low, in some cases involving transition and relaminarization. In 
these cases, the thickness of the laminar sublayer becomes an appreciable fraction of the 
overall boundary layer thickness and the wall function approach becomes inapplicable or 
ineffective. In addition, the wall function approach gives questionable results for separated 
flows. For these reasons it was decided to use the wall damping function approach. 

Two equation models have been emphasized since these are viewed as the sim- 
plest and most practical models available which have sufficient generality to be applied to 
the complex flows of interest in NASP applications. Although numerous two equation 
models have been investigated, descriptions of only two of these will be given here since 
these are the baseline models of our final recommendation. In Coakley and Huang (1992) 
a detailed investigation of the performance of a variety two equation models was pre- 
sented for flat plate boundary layers over a wide range of Mach and Reynolds numbers. 
The results of that study showed that most of the models gave reasonably good predictions 
of skin friction, heat transfer and velocity profiles with little clear preference of one model 
over the other. For this reason only two models were selected for further study, and it is 
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believed that these models are representative of most two equation models currently in 
use. The baseline models are the k-e model of Jones and Launder (1972), as modified by 
Launder and Sharma (1973), and the k-at model of Wilcox (1984). These models are used 
with the (mass weighted) Reynolds averaged compressible Navier-Stokes equations 
which, in cartesian tensor form, are given below. 

Reynolds Averaged Compressible Navier-Stokes Equations 

|(p) + 4 <p " ; ' ) =0 

^ (p "' ) + ^. <p u ‘ u i + °i? =0 

& (p£)+ 5f/ p =° 


where; p is the density; « t - are the cartesian velocity components; E = e + 05u{Ui + k is the 
total specific energy; e = c v T is the specific internal energy; k-05 pufuf /p is the turbu- 
lent kinetic energy; T is the temperature; p = (y - /) pe is the equation of state; and p is 
the pressure; y = Cp/c v is the ratio of specific heats, and c p and c v are the specific heats at 
constant pressure and volume respectively. The variables Oj- and qj are the total stress 
tensor and heat flux vector, respectively, which include both molecular and (Reynolds 
averaged) turbulent contributions. Using the Boussinesq approximation, these variables 
are represented in terms of an eddy viscosity by 
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where H and jJ. r are the molecular and turbulent (eddy) viscosities a and oj are molecular 
and turbulent Prandtl numbers with a = c p \l/K (assuming air) a T = 0.9 and ct* depending 
on the model used. The turbulent eddy viscosity is expressed in terms of the turbulent 
kinetic energy, k, and either the dissipation rate, e, or the specific dissipation rate o> 
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depending on the model. This expression is 


Mr - Cyf^pjlcl - p Cyfp P 



where Vit is the turbulent velocity scale, / = W/£ = Va/cd is the length scale, is a mod- 
eling constant, and is a damping function depending on the specific model used. For 
applications to complex flows, the governing equations were expressed in terms of curvi- 
linear coordinates and solved using the finite volume method, Viegas and Rubesin (1991), 
and Huang and Coakley (1992b) or the finite difference method, Bardina (1994). 

The recommended baseline k-e and k - co turbulence models are expressed by the 
formulas given below. 


3(p*, + ^(p ku r » t *) = port 
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In these equations and in the following Tables, variable “s” and subscript “s” are replaced 
with e for the k-e model or co for the k-(o model, respectively. 

The model parameters are defined in Tables 1, 2 and 3 


Table 1: Model Parameters 


II 
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Table 2: Launder-Sharma k-e Model, s = e = cak 


o k = 1 

a £- 1-3 

cP 

II 

u\ ! 

C e2 = 1.92 

= exp( - 3.4/(1.0 + R 2 /50)) 

h = 1 - 0.3 exp( - R 2 ) 

D k =l + (2v/e) (dk ll2 /dx k ) 2 

d e = Q2/2 - ( 2v Ve 2 ) (^(u^^/dxfdx ^ 2 


Table 3: Wilcox Jt-co Model, s ■ co = efk 


o k = 1/2 

c^co = 1/2 

c a , = 519 

c (a2 = 5/6 

H 

II 

II 

>— * 

D a - C(o2 


For the Launder-Sharma model, the thin layer approximation was normally used 
to compute the derivatives in the D k and D e terms. 

An important consideration in using the Wilcox k - co model, which is not neces- 
sary with the k-e model, is the value of co in the free stream (just outside the boundary- 
layer edge) where it cannot be too small. In all of the applications of this study, it was pos- 
sible to choose values of co at the inflow boundary which insured that the values of co in 
the free stream would not be too small. (It must be chosen such that tOoo > 10 U<JL where 
Uoo is the free stream velocity and L is the length of run of the boundary layer and (Ooo is 
the free stream value of co at the start of the boundary layer). It may be that in future appli- 
cations it will not be feasable to control the free stream C0oo in this manner, and other mea- 
sures will be necessary. One alternative would be to use the k - co model of F. Menter 
(1992), which uses a blending of the Jt-co and k-e models to circumvent the problem. This 
model will be discussed more fully in the next section. 

Recommended Model Corrections for Complex Flows 

The baseline models described above are generally not adequate to accurately 
predict complex flow problems and must be corrected to deal with these cases. The types 
of modeling corrections that have been found useful in practice and which are recom- 
mended for NASP applications are summarized below. Other model corrections which 
have been tried in the course of the study but which are not recommended for applications 
are described in Coakley and Huang (1992). 
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Length Scale Correction 


The first correction is addressed to difficulties encountered in predicting heat 
transfer in the reattachment or shock impingement zone of shock-wave boundary-layer 
interaction flows. In these zones, all 2-eq models dramatically overpredict heat transfer 
and must be corrected. The correction involves the use of an algebraic length scale which 
limits the slope of the length scale predicted by the two equation model, which otherwise 
would become very large in these regions. The formulas defining the length scale correc- 
tion for the k-e and k-(D models are given below, 

l = min{2.5y,Jk?/E] = min{2.5y,Jk/a>} 

in these formulas, / is the turbulent length scale which is taken to be the smaller of an alge- 
braic expression (k y — 2.5 y, based on a von Karman constant of k = 0.41) and the 
conventional length scale given by the two equation model. Having computed this param- 
eter, the value of e or cd is recomputed and reset to be consistent with this value, e.g. e = 

V* 3 // or co= V*//. 


Rapid Compression Correction 

The second correction to be described is called the rapid compression correction 
and is used to improve predictions of separation in shock-wave turbulent boundary-layer 
interactions. The correction involves changing the coefficient of the dilitation or velocity 
divergence in the e and co equations, i.e. or a w The net effect of this correction is to 

increase the production of epsilon or omega in regions of rapid compression, or shock in 
waves, which reduces the eddy viscosity and enhances separation. The corrected values of 
the dilitation coefficient for each model are 

a e = 2 or = 4/3 

as opposed to the Launder-Sharma coefficient a E = (2/3) Q/ = 0.97 and the Wilcox coeffi- 
cient a w = (2/3) C mI = 0.37 shown in Tables 1, 2, and 3. The development of this correc- 
tion is discussed more fully in Coakley and Huang (1992). 
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Compressible Dissipation Correction 

The third and final correction recommended for NASP applications is applicable 
to free shear flows (e.g. mixing layers, wakes and jets). In these types of flows, it is well 
known that shear layer spreading rates decrease as Mach numbers increase compared with 
spreading rates at zero Mach number. The correction recommended to improve predic- 
tions of these flows was developed originally by Zeman (1990), and is closely related to 
similar corrections developed by Sarkar (1991), and Wilcox (1992). The correction is 
listed below, 

D lc D k + 

D CD D co ” ^ 

<j) fc = a 0 (l - exp (-(max (0, a^j- a 2 ) /a 3 ) 2 )) 

D e is unchanged, M r - 'Jk/c, c is the local sound speed, ag = 3/4, aj = V(y-t-l), a 2 = 1/10, 
and a$ - 6/10. 

For applications involving boundary layers, or solid walls, this correction has 
been found to underpredict skin friction especially at high free stream Mach numbers. It is 
recommended, therefore, that this correction only be used in free shear flow applications. 
The correction of Wilcox (1992), which is similar to that of Zeman, has been found to 
work well in boundary layers as well as free shear flows and Wilcox recommends its use 
without reservation. Since we have not investigated this model under the wide range of 
conditions investigated using other model corrections, we chose not to recommend the 
model at this time. 


Other Corrections and Models of Interest 

Although the above models and corrections constitute our final recommenda- 
tions for NASP applications, it is important to note that these models only constitute an 
improvement over previous models and may very well give poor predictions for flow situ- 
ations which have not yet been investigated and validated. In this regard we draw attention 
to the k-CD model of F. Menter (1993), which has proven quite successful in incompress- 
ible and transonic flow applications. This model was designed to overcome certain defi- 
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ciencies of the Wilcox k-co model and has certain features which might prove useful in 
hypersonic applications. One feature of the model is the so called shear stress transport or 
rapid strain feature which enables the model to give improved predictions of adverse pres- 
sure gradient boundary-layers and separation. This feature is similar to the rapid compres- 
sion correction discussed above and enhances separation by reducing the eddy viscosity in 
non-equilibrium regions where the flow is changing rapidly. This feature was also tried 
with the baseline models and gave results very similar to those obtained with the rapid 
compression correction. It is believed that some combination of the two corrections may 
ultimately prove more accurate and reliable in future applications. It must be stated, how- 
ever, that we did apply Menter’s model to most of the flows described in this study and 
obtained no improvement over the recommended models. In some cases it did not perform 
as well. Since the model is considerably more complicated than the other models studied 
we decided not to include it in the list of recommended models. 

Model Validation - Representative Results 

Representative results of calculations and comparisons of model predictions with 
experimental measurements will be given in this section. The flows discussed include free 
shear (mixing layer) flows, 3 two dimensional shock- wave boundary- layer interaction 
flows, and 2 three dimensional shock-wave boundary-layer interactions. The turbulence 
models used include the recommended baseline and corrected k-£ and k-co models. 

Compressible Mixing layer 

The first series of flows to be discussed consists of high speed mixing layers 
which are of considerable importance in the design of propulsive exhaust nozzles for the 
NASP vehicle. The comparisons of computations with experimental measurements is 
shown in fig.(5). The figure shows predictions of spreading rate divided by spreading rate 
at zero Mach number compared with experimental measurements over a range of convec- 
tive Mach numbers. The data include the Bogdanoff (1993) compilation of the Langley 
data, and the measurements of Samimy and Elliot (1990). The calculations were done by 
Viegas and Rubesin (1992), who used the baseline k-e model with the compressible dissi- 
pation corrections of Zeman and Sarkar. Although not shown, results obtained with the k- 
© model give similar results. Examination of the comparisons indicates that the baseline 
model significantly overpredicts the spreading rate while the model corrections improve 
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the predictions. Of the two model corrections, the Zeman correction gives the best overall 
result. 

Two-Dimensional Shock -Wave Boundary-Layer Interaction Flows 

The second series of model validation studies to be discussed consists of two 
dimensional (planar and axisymmetric) shock- wave boundary-layer interactions. The 
flows are the Mach 7 ogive-cylinder- flare flow of Kussoy and Horstman (1989) (35° flare 
angle), the Mach 9 planar compression ramp flow of Coleman and Stollery (1972) (34° 
ramp angle), and the Mach 7 axisymmetric impinging shock flow of Kussoy and Horst- 
man (1975) (15° generator angle). In all cases, the walls were highly cooled with wall-to- 
adiabatic wall temperature ratios on the order of 0.3 to 0.4. The test configurations for 
these cases are shown in figs.(6a,b,c). 

All calculations were done using the code developed by Huang and Coakley 
(1992). The inlet flow conditions just ahead of the shock interaction zones were obtained 
by calculating the flow over a flat plate and matching measured and computed displace- 
ment thicknesses. The value of y + at the first grid point off the wall was maintained to be 
less than 0.5 and the grid was expanded exponentially from the wall to the free stream. 
This gave between 60 to 80 grid cells in the boundary-layer and 140 cells overall in the 
cross stream direction. Computations with fewer cells inside the boundary-layer (i.e. 40 
cells) were made and no significant differences were observed. In the streamwise direction 
a grid of 140 cells was used except in the impinging shock case where 200 cells were used. 
The comparisons of computations with measurements are shown in figs.(7-9). They 
include measured and computed surface pressures and heat transfer distributions (and skin 
friction for the impinging shock case) and were done with the baseline and corrected k-e 
and k-<£> models. The pressure and heat transfer measurements are normalized by the mea- 
sured values in the region ahead of the interaction. Zeman ’s compressible dissipation cor- 
rection designed for shear layers was not used as explained in the section on model 
corrections. It is clear from these results that the baseline models significantly under pre- 
dict the extent of separation and over predict the heat transfer in the interaction region. 
The corrected models both give results in much better agreement with experiment. 
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Three-Dimensional Shock-Wave Boundary-Layer Interaction flows 

The Ames experiments on 3-D shock-wave boundary-layer interactions 
described in this report were used here to test the compressible turbulence models and the 
model corrections. The experiments and the data are also described in Settles and Dodson 
[1993a and 1993b] and Kussoy and Horstman [1991,1992,1993]. The calculations were 
done with the code of Bardina (1994) 

3-D Vertical Fin Shock Interaction Flows 

This experiment investigates the interaction of a hypersonic shock wave with a 
thick turbulent boundary layer [Kussoy and Horstman, 1991 and 1993b]. A 10° and a 15° 
vertical fins mounted on top of a flat plate were used to generate oblique shock waves. The 
free-stream Mach number was Moo = 8.2, the temperature was Too = 81° K, and the Rey- 
nolds number was Reoo = 5- 10 6 per meter. The wall temperature was fixed at 300° K. The 
interaction of the shock wave with the turbulent boundary layer generates a crossflow vor- 
tex separation with a “quasi conical” shape [Settles and Lu, 1985; Knight, Horstman, and 
Monson, 1992], Peak wall pressure, skin friction, and heat transfer rates were observed in 
the re-attachment zone behind the crossflow vortex. 

The numerical computations were made with 61x41x61 an 31x21x31 meshes 
[Bardina, Coakley, and Marvin, 1992]. Only small differences between the solutions were 
observed, and the fine mesh solutions are considered accurate for engineering purposes. 
The inflow conditions were obtained from the Navier-Stokes code solution matching the 
experimental displacement thickness. 

A few comparisons of experiment and simulation with the finer mesh are 
described below. Figures 10a and 10b show the surface pressure and skin friction distribu- 
tions, respectively, for the 10° fin flow on the flat plate surface at the crossed section 
located at x=0.1819 m downstream of the fin leading edge. Figure 10c shows the compar- 
ison of the wall heat transfer distribution on the flat surface at x=0.1645 m downstream of 
the fin leading edge. Comparable results for the 15° fin case are shown in fig.(lla,b,c). 
The symbols in the figures show the experimental data points, the solid lines show the 
solution with the baseline k-a model, and the dash lines show the solution with the k - co 
model with both model corrections (length-scale and rapid compression corrections). Both 
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simulations show good agreement with the experimental data. The peak values observed 
in the re-attachment zone are also well predicted. The simulations fail to predict the small 
plateau observed in the wall heat transfer rate distributions upstream of the shock wave. 
Comparison of baseline and corrected model predictions show only small differences. 
This is probably due to the fact that, compared with the 2-D results, the pressure rise 
through the shock wave is relatively weak and separation is relieved by three dimensional 
effects. 


3-D Crossing Shock Interaction Flow 

This experiment studies the interactions of two intersecting hypersonic shock 
waves with a thick turbulent boundary layer [Kussoy and Horstman, 1992], Two 15° fins 
mounted on top of a flat plate were used to generate intersecting oblique shock waves (see 
fig. 4). The free-stream Mach number was Me = 8.3, the temperature was T TO = 80° K, and 
the Reynolds number was Reoo = 5.3' 10** per meter. The wall temperature was fixed at 
300° K. The intersection of both crossflow vortices generates different complex flow 
structures with high static pressures and surface heat transfer rates. The intersection of the 
two “quasi conical” vortical structures uplifts the flow and induces a wave structure in the 
symmetry plane [Gaitonde and Shang, 1993]. 

The numerical simulations were made with a 231x81x81 mesh [Bardina and 
Coakley, 1994; Bardina, 1994], Simulations studies with 101x61x41 and 31x21x31 grid 
points were also done to analyze grid effects. Small differences between the solutions 
were observed in the surface pressure and heat transfer distributions, but significant dif- 
ferences were observed in flow structure. The fine mesh solutions provided the best reso- 
lution of the turbulence structures and are considered accurate for engineering purposes. 
The inflow conditions were obtained from the Navier-Stokes code solution matching the 
experimental displacement thickness. 

Figures 12a and 12b show the pressure and heat transfer distributions, respec- 
tively, on the plate surface along the symmetry plane located between the two fins. The 
symbols show the experimental data, the solid lines show the solution with the baseline k- 
Q) model, and the dash lines show the solution with the k < o model with both model correc- 
tions (length-scale and rapid compression corrections). Predictions of both surface pres- 
sure and heat transfer show good agreement, within the experimental uncertainty, except 
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near the outflow zone. The small plateau at the beginning of the interaction is also not pre- 
dicted by either model. The peak pressure and heat transfer rate are very well predicted. 
As in the single fin cases, both the baseline and corrected model predictions show only 
small differences. 

Figures 13a and 13b show the pressure distributions on the plate surface at two 
cross-sections, one located at x/Sm = 5.60 and the other located at x/S,*, = 6.92 downstream 
of the fin leading edge. The first distribution is upstream of the peak surface pressure, and 
the second distribution is downstream of the peak surface pressure generated by the flow 
re-attachment. Both simulations show good agreement within the experimental uncer- 
tainty of the data, and both models give similar predictions. 

Figures 13c and 13d show the comparison of heat transfer rate profiles on the 
plate surface in the crossed sections located at x/5oo — 5.08 and x/5,* = 6.40, respectively, 
upstream and downstream of the peak heat transfer rate observed in the symmetry plane. 
Both simulations agree in general with the experimental data. The model corrections show 
improved heat transfer rate predictions in the downstream zone. 

Figure 14 shows experimental pitot pressure contours compared with computa- 
tional contours obtained with the corrected k-co model. Three locations are shown, = 
5.60, 6.9, and 8.3, respectively. The resolution of the numerical data is 81x81 grid points 
while the resolution of the experimental data is 5x24, 4x24, and 4x24 respectively. The 
agreement with the experimental data is very good. 

Summary and Conclusions 

In this section we summarize the work performed under GWP 18, give our prin- 
cipal results and recommendations, and discuss plans for future work. We feel that, over- 
all, the results produced in the course of the work were of a very high caliber and will be 
of considerable use to modelers and designers of hypersonic flight vehicles. The research 
was divided into two distinct phases; one experimental and the other theoretical. The 
experimental phase consisted of the collection and assessment of a database of high speed 
wind tunnel experiments gathered from sources around the world and the conduct of addi- 
tional experiments in the NASA Ames 3.5’ hypersonic wind tunnel which would add to 
the database. The primary aim of this effort was to produce a database of reliable and rele- 
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vant flow measurements which could be used to validate turbulence models for NASP 
CFD design codes. Out of a total of over 2000 experiments that were initially studied, the 
list of acceptable experiments was reduced a total of 30 using a set of strict acceptance cri- 
teria. The kinds of flows analyzed and recommended as legitimate candidates for inclusion 
in the database included two- and three-dimensional shock-wave boundary-layer interac- 
tion flows, attached boundary layers flows and free shear flows. 

The theoretical phase of the work consisted of identifying, testing, and recom- 
mending turbulence models which would be of practical use in the CFD design codes. A 
large number of baseline models were initially tested on flat plate flows and this number 
was then reduced to two for further testing on more complex flows. The baseline models 
selected were the Launder-Sharma version of the k-e model, and the Wilcox k - co model. 
These models were tested on free shear flows and two- and three-dimensional shock-wave 
boundary-layer interaction flows. Some of the results of these calculations have been dis- 
cussed in this report. It was found that the baseline models did not perform satisfactorily 
with regard to separation and heat transfer predictions especially for the 2-D shock-wave 
boundary-layer interaction flows, and did not accurately predict the spreading rate of free 
shear layers. To improve model performance for these complex flows, a series of com- 
pressibility corrections was investigated. The more promising of these corrections were 
then selected to be the final model recommendations for incorporation in to the NASP 
CFD codes. Representative results using the corrected models were discussed in this 
report. It was shown that for the flows investigated, the model corrections give substan- 
tially improved predictions. 

Although the work performed in the course of this research has led to the identi- 
fication of useful flows for model validation and the development of improved turbulence 
models for hypersonic flight, much work remains to be done. Not all of the flows included 
in the database of recommended flows have been investigated computationally, and these 
need to be investigated. More complicated flows not included in the database also need to 
be investigated. Flows that fall into this latter category include flows with chemical reac- 
tions such as those occurring in SCRAM jet combustion and propulsion, and flows that 
include transitional phenomena just to mention two. Interest in these and other areas of 
turbulence modeling research is very high and the work goes on. 
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Fig. 5. Effect of compressibility on normalized vorticity thickness growth rate for mix- 
ing layer (Viegas and Rubesin (1991, 1992)). 


27 



M = 7.0 



Fig. 6a. Axisymmetric compression generated with ogive-cylinder-flare configuration of 
Kussoy and Horstman (1989). 
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Fig. 6b. 2-D planar compression ramp configuration of Coleman and Stollery 

(1972). 
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Fig. 6c. Axisymmetric impinging shock generated with a cone-ogive cylinder and an 
annular fin configuration of Kussoy and Horstman (1975). 
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Fig. 7. Surface pressure (a) and heat transfer (b) in axisymmetric compression 
flow of Kussoy and Horstman (1989), geometry shown in Fig. 6a. 
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Fig. 8. Surface pressure (a) and heat transfer (b) in 2-D planar compression ramp 
flow of Coleman and Stollery (1972), geometry shown in Fig 6b. 
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Fig. 10. Surface flow predictions in the 3-D vertical fin interaction flow (a=10° 
case) of Kussoy and Horstman (1991a), geometry shown in Fig. 3. 
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Fig. 1 1 . Surface flow predictions in the 3-D vertical fin interaction flow (a=15° 
case) of Kussoy and Horstman (1991a), geometry shown in Fig. 3. 
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Fig. 12. Surface flow predictions in the symmetry plane of the 3-D crossing shock 
flow (ct=15° case) of Kussoy and Horstman (1992), geometry shown in Fig. 4. (Symbols 
represent experimental data, — k - © model, and — k-ob model with corrections). 
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Fig. 13. Surface flow quantities in the 3-D crossing shock interaction flow (a=15 
case) of Kussoy and Hors tm an (1992), geometry shown in Fig. 4. (Symbols represent 
experimental data, — k - <d model, and — k - © model with corrections). 
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Fig. 14. Pitot pressure contours on cross-sections of 3-D crossing shock interaction 
flow (cc=15° case) of Kussoy and Horstman (1992), geometry shown in Fig. 4. (The fin 
leading edges are located at x = 0 , and 800 is the thickness of the undisturbed turbulent 
boundary layer at x = 0 ). 
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